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Abstract 

Synchronization of coupled harmonic oscillators is investigated. Cou- 
pling considered here is pairwise, unidirectional, and described by a non- 
linear function (whose graph resides in the first and third quadrants) of 
some projection of the relative distance (between the states of the pair 
being coupled) vector. Under the assumption that the interconnection 
topology defines a connected graph, it is shown that the synchronization 
manifold is semigiobally practically asymptotically stable in the frequency 
of oscillations. 



1 Introduction 

Synchronization in coupled dynamical systems, due to the broad range of appli- 
cations, has been a common ground of investigation for researchers from different 
disciplines. Most of the work in the area studies the case where the interconnec- 
tion between individual systems is linear; see, for instance, [TT1 151 [TT1 [Tbl HI 17] . 
Nonlinear coupling is also of interest since certain phenomena cannot be mean- 
ingfully modelled by linear coupling. A particular system exemplifying non- 
linear coupling that attracted much attention is Kuramoto model and its like 
[5J [TD]. Among more general results allowing nonlinear coupling are [HQS] where 
passivity theory is employed to obtain sufficient conditions for synchronization 
under certain symmetry or balancedness assumptions on the graph describing 
the interconnection topology. 

In this paper we consider the basic equation of the theory of oscillations 

q = -~uj 2 q 

i.e., the harmonic oscillator. One physical example is a unit mass attached to a 
spring. We interest ourselves with the following question. If we take a number of 
identical mass-spring systems and couple some pairs by unidirectional nonlinear 
dampers, will they eventually oscillate synchronously? To be precise, do the 
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solutions of the following array of coupled harmonic oscillators 
q\ = -uj 2 q t + ^2 7w i = 1) 2 ) • • • , P 

(where jij(-) is either identically zero or it is some function whose graph lies 
entirely in the first and third quadrants) synchronize? This question may be of 
direct importance for certain simple mechanical systems or electrical circuits, 
but our main interest in it is due to the possibility that it may serve as a lucky 
starting point for understanding a more general scenario. 

Regarding the above question, our finding in the paper is roughly that syn- 
chronization occurs among mass-spring systems if the springs are stiff enough 
and there is at least one system that directly or indirectly effects all others. More 
formally, what we show is that, for a given set of functions {7y(-)} describing the 
coupling configuration, if the graph representing the interconnection topology 
is connected, then the solutions can be made to converge to an arbitrarily small 
neighborhood of the synchronization manifold, starting from initial conditions 
arbitrarily far from it by choosing large enough u>. In technical terms, what 
we establish is the semiglobal practical asymptotic stability (in uj) of the syn- 
chronization manifold. Intuition and simulations tell us that global asymptotic 
synchronization should occur regardless of what to is. This however we have not 
been able to prove (nor disprove). 

We reach our final result in three steps. Note that the solution of an un- 
coupled harmonic oscillator defines a rotating vector on the plane. Thanks to 
linearity of the system the speed of this rotation (ui) is independent of the ini- 
tial conditions. As a first step therefore we express the systems with respect to 
a rotating coordinate system. This change of variables yields coupled systems 
whose righthand sides are periodic in time (with period 2tt/lu). Our second step 
is to exploit this periodicity. We obtain the average systems and realize that 
they belong to a well-studied class of systems pertaining to consensus problems 
[5] . We then deduce that the solutions of average systems converge to a fixed 
point on the plane. As our final step we use the result of [H] to conclude that 
the global asymptotic synchronization of average systems implies the semiglobal 
practical asymptotic synchronization of coupled harmonic oscillators. 

2 Preliminaries 

Let N denote the set of nonnegative integers and R>o the set of nonnegative 
real numbers. Let | ■ | denote Euclidean norm. For 1 with 

Xi 6 R™ we let A := {x £ M. np : xt — Xj for all i, j} be synchronization 
manifold. A function a : K>o — > K>o is said to belong to class-/C (a 6 fC) if it 
is continuous, zero at zero, and strictly increasing. A function /3 : K>o X K>o — > 
M>o is said to belong to class- ICC if, for each t > 0, (3(-,t) is nondecreasing 
and lim s ^ + [3{s, t) — 0, and, for each s > 0, (3{s, ■) is nonincreasing and 
limt-xxj (3(s, t) = 0. Given a closed set S C l n and a point x £ R™, |a;|s 
denotes the (Euclidean) distance from x to S. 
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A (directed) graph is a pair (j\f, £) where Af is a nonempty finite set (of nodes) 
and £ is a finite collection of ordered pairs (edges) (rii, rij) with rii, rij G Af. A 
directed path from ni to is a sequence of nodes (m, n%, ■ ■ • , n<) such that 
(m, rij+i) is an edge for i G {1, 2, ...,£— 1}. A graph is connected if it has a 
node to which there exists a directed path from every other node0 

A set of functions {7.^ : M — > R}, where i, j' = 1, 2, . . . , p with i 7^ j, 
describes (is) an interconnection if the following hold for all i, j and all s£l: 

(i) 7ij (0) = and 37^ (s) > 0. 

(ii) Either 7y(s) = or there exists a G K. such that |7y(s)| > a(|s|). 

To mean "fij(s) = we write 7^ = 0. Otherwise we write 7^ ^ 0. The graph 
of interconnection {7^} is pair (A/", £), where Af = {ni, . . . , n p } and £ is such 
that (ni, nj) £ £ iff 7^ ^ 0. An interconnection is said to be connected when 
its graph is connected. 

To give an example, consider a set of functions C := {7^ : i, j = 1, . . . , 4}. 
Let 713, 723, 724, 732 be as in Fig. 1 while the remaining functions be zero. 
Note that each 7^ satisfies conditions (i) and (ii). Therefore set C describes an 
interconnection. To determine whether C is connected or not we examine its 
graph, see Fig. 2. Since there exists a path to node 77,4 from every other node, 
we deduce that the graph (hence interconnection C) is connected. 

3 Problem statement 

We consider the following array of coupled harmonic oscillators 



where uj > and {7^ } is a connected interconnection. We assume throughout 
the paper that 7^ are locally Lipschitz. Let £j G US 2 denote the state of ith os- 
cillator, i.e., £j = [qi Pi] T ■ When coupling functions 7^ are linear, the oscillators 
are known to (exponentially) synchronize for all uj. That is, solutions £i(-) con- 
verge to a common (bounded) trajectory, see [T5]. In this paper we investigate 
the behaviour of oscillators under nonlinear coupling. In particular, we aim to 
understand the effect of the frequency of oscillations uj on synchronization for a 
given interconnection {jij}. 





(la) 
(lb) 



4 Change of coordinates 



We define S(uj) G R 2x2 and H G R lx2 as 



uj 
-uj 



H := [0 1] . 



This is another way of saying that the graph contains a spanning tree. 
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Figure 1: Some examples of coupling functions. 



Then we rewrite (JXJ) as 

it = S(w)ti + H T £ 7 «(ff & - SO) ■ (2) 

Let us recall the geometric meanings of the terms in @. The first term S(u>)£i 
defines a rotation (with period 27r/w) since S(u>) is a skew-symmetric matrix. 
The second term H T — &)) induces a contraction in the direction 

specified by vector H T (so to speak, along the vertical axis) such that the 
projections of states £j on the vertical axis tend to approach each other. The 
combined effect of these two terms is relatively more difficult to visualize. One 
trick to partially overcome this difficulty is to look at the system from the point 
of view of the observer that sits on a rotating frame of reference. 

Under change of coordinates Xi(t) := e we can by J2J write 

V^7ij([~ sin uit cos u)t](xj — xi)) . (3) 
Since the change of coordinates is realized via rotation matrix e~ s ^ 1 , the 



— sin tut 
cos ut 
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Figure 2: Graph of interconnection C. 



relative distances are preserved, that is \xi(t) — Xj(t)\ — |£i(t) — £j(t)\ for all t and 
all i, j. This means from the synchronization point of view that the behaviour 
of array @ will be inherited by array ©. However, exact analysis of seems 
still far from yielding. Therefore we attempt to understand this system via its 
approximation. 



5 Average systems 



tions 7ij : R 2 



Observe that the righthand side of ([3]) is periodic in time. Time average func- 

7ij([— shop cos ip]x)d(p . 



7«(*) ■■= ± 



;iven 


by 








— sin ip 


Jo 


COS if 



Then the average array dynamics read 

m = J2iij(vj -Vi)- ( 4 ) 

Theory of perturbations [3j Ch. 4 § 17] tells us that, starting from close initial 
conditions, the solution of a system with a periodic righthand side and the solu- 
tion of the time-average approximate system stay close for a long time provided 
that the period is small enough. Therefore ((4]) should tell us a great deal about 
the behaviour of ^ when uo ^> 1. 

Understanding requires understanding average function 7^. The follow- 
ing lemma is helpful from that respect. 

x 

Lemma 1 We have Jij(x) = Pij(\x\)-. — r where p t j : R> — > R>o *s 

\x\ 

1 f 2 * 

Pij( r ) '■= 7T" / Jij {r sin (p) simp d<p. 
27T Jo 

Proof. Given x e R 2 , let r = \x\ and € [0, 2n) be such that r[— cos 8 sm8] T = 



5 



x. Then, by using standard trigonometric identities, 

7*j (''(sin tp cos 9 + cos tp sin 6)) dip 
7ij(rsin(</5 + 0))dy> 



Z7T 



- sin ip 
cos 







1 

2^ 
1 

2^ 
1 

2Wo 
1 



- sm (/? 
cosy; 



2n 



— s'm((p — 9) 
cos(ip — 9) 



sm tp cos ft + cos tp sm 
cos </? cos 9 + sin sin 



7ij(r sintp)^ 

7y (rsint^)^ 



271 



= ( 2~ / 7jj( rsin ^)sin(y5d^ 



cos ( 
sin( 



1 

2^ 
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7y (r sin tp) cos </? dtp 



sin (9 
cos# 



(5) 



We focus on the second term in ([5]). Observe that 

f-n/2 



f jij (r sin tp) cos tp dtp = J ji 3 ■ (r sin ^ + cos ^ + ^ d<y3 
= 

since the integrand is an odd function on the interval of integration. Likewise, 



jij (r sin tp) cos 



tt/2 



-tt/2 



3tt 



7y r sin [tp + — cos Uo + — ) dtp 



3- 



= 0. 



Therefore 



jij (r sin tp) cos tp dtp 



jij (r sin tp) cos tpdtp + / 7^ (r sin tp) cos </j 



0. 



(6) 



Combining ([5]) and ([5]) we obtain 



Jij(x) = 



1 

2^ 



7y (r sin sin tp dtp 



— cos ( 
sin( 



= Pij(\x\h 



Hence the result. 
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Lemma [T] tells us that, given any vector x on the plane, vector 7y(a:) (if 
it has nonzero magnitude) is in the same direction as x. In the light of this 
if we now look at systems (0| we can roughly visualize the evolution of the 
trajectories. Take the ith system. Let ji, . . . , jt be all the indices such that 
j G {ji , . . • , je} implies 7^ ^ 0. These indices are sometimes called the indices 
of neighbors of system i. Then we can write f)i = Vi + . . . + vt where each Vk is a 
vector pointing from rji to the state of fcth neighbor system. Therefore the net 
velocity vector ^ Vk points to some "weighted mean" of the neighbor systems' 
states. Synchronization of such systems, where the velocity vector of a system 
always points to some weighted mean of the positions of its neighbors, have 
been studied under the names consensus and state agreement] see, for instance, 
[H [TJ [6] . The finding of those works is roughly that if the interconnection is 
connected, then the solutions of systems converge to a common fixed point in 
space. We now give the formal application of this result to our case. 



Let us stack the individual states r\i to form r\ := [rj[ r\^ 



TF. Define 



7av(?7) 

We then reexpress (UJ) as 



E7ii(%- -m) 

V = 7avfa) • 



(7) 



We now have the following result. 



Theorem 1 Consider system ([7|). Synchronization manifold A is globally asymp- 
totically stable, i.e., there exists /3 G K-C such that \t](t)\^ < /3(|^(0)|^, t) for all 
t>0. 



Proof. By Lemma Q] we can write 



m\ 



for % = 1, 2, . . . , p. Let Q denote the graph of interconnection {7^ }. Note 
that Q is connected by assumption. We make the following simple observations. 
Function pij is continuous and zero at zero. If there is no edge of Q from node 
i to node j then Pij{r) = 0. If there is an edge from node i to node j then 
(7jj (r) > for r > 0. 

Therefore fi is continuous; and vector fi(rj) always points to the (relative) 
interior of the convex hull of the set {?7i}U{rjj : there is an edge of Q from node i 
to node j}. These two conditions together with connectedness of Q yield by 
[51 Corollary 3.9] that system ([7]) has the globally asymptotic state agreement 
property, see [6j Definition 3.4]. Another property of the system is invariance 
with respect to translations. That is, 7 av (?7 + = 7av(?7) f° r £ £ A. These 
properties let us write the following. 
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(a) There exists a class-/C function a such that \r)(t)\j( < a(|r/(0)|^) for all 
t > 0. 

(b) For each r > and e > 0, there exists T > such that |t^(0)|^ < r implies 
\v(t)U < e for a11 f > T - 

Finally, (a) and (b) give us the result by [T51 Proposition 1]. ■ 

Let us now go back to our discussion in the beginning of Section [4] There 
we talked about two actions that shape the dynamics of systems @, namely, 
rotation and vertical contraction. The combined effect of those actions on syn- 
chronization of the systems was not initially apparent. However, by applying 
first a change of coordinates and then averaging (HJ) we see that it is likely 
that two actions will result in synchronization, at least when the rotation is 
rapid enough. Vaguely speaking, rotation rescues contraction from being con- 
fined only to vertical direction and sort of smears it uniformly to all directions, 
which should bring synchronization. In the next section we formalize our ob- 
servation. 



6 Semiglobal practical asymptotic synchroniza- 
tion 



Consider systems ([3]). Stack states Xi to form x 



„TlT 



Define 



7(x, Lot) := 



- sin uit 
cosojt 

- sin cot 

COS LOt 



S 7ij([~ siniot cos u>t](xj — xi)) 



j p j([— sin Lot cosu>t](xj — x p )) 



Now reexpress ^ as 

x = 7(x, LOt) . 

The following definition is borrowed with slight modification from [14] . 



(8) 



Definition 1 Consider system x = /(x, Lot) . Closed set S is said to be semiglob- 
ally practically asymptotically stable if for each pair (A, S) of positive numbers, 
there exists lo* > such that for each lo > lo* the following hold. 

(a) For each r > S there exists e > such that 

|x(0)| s <£ |x(t)| s <r Vi >0. 

(b) For each e < A there exists r > such that 

|x(0)| 5 <e |x(t)| s <r Vi >0. 
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(c) For each r < A and e > 6 there exists T > such that 



|x(0)| s <r 



|x(t)| s <e Vt>T. 



Below we establish the semiglobal practical asymptotic stability of synchroniza- 
tion manifold. To do that we use [TH Thm. 2], which says that the origin of 
system x = /(x, cot) (where / is periodic in time) is semiglobally practically 
asymptotically stable if the origin of x = / av (x) (where / av is the time average 
of /) is globally asymptotically stable. 

Theorem 2 Consider system (|5J|. Synchronization manifold A is semiglobally 
practically asymptotically stable. 

Proof. Consider systems (0) . Observe that the righthand side depends only on 
the relative distances Xj — Xj. Let us define 



x 2 - Xi 
x 3 - Xi 

y := 

- xi 

Note that y = f(y, ujt) for some / : M 2 ^ 2 x R> -> R 2p " 2 . Since functions 
jij are assumed to be locally Lipschitz, / is locally Lipschitz in y uniformly 
in t. Also, / is periodic in time by ((3]). Now consider systems (J4j) . Again the 
righthand side depends only on the relative distances rjj — Define 

V2 - m 
m - m 

n P - m 

Then z = / av (z) where / av is the time average of / and locally Lipschitz both 
due to that 7^ is the time average of 7^ . 

Theorem Q] implies that the origin of z = f av (z) is globally asymptotically 
stable. Then [TH Thm. 2] tells us that the origin of y = f(y, cut) is semiglobally 
practically asymptotically stable. All there is left to complete the proof is to real- 
ize that semiglobal practical asymptotic stability of the origin of y = /(y, cot) is 
equivalent to semiglobal practical asymptotic stability of synchronization man- 
ifold A of system ■ 

Theorem [5] can be recast into the following form. 

Theorem 3 Consider coupled harmonic oscillators (fTJ). For each pair (A, 5) 
of positive numbers, there exists co* > such that for each co > co* the following 
hold. 
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(a) For each r > S there exists e > such that 

max |6(0) - £,*(0)| < £ => max |6(t) - £,•(*) | < r Vt > . 



(7>) For each e < A £/iere exists r > swc/i t/wrf 

max |6(0) - 6/(°)l < £ =^ max - &(*)| < r Vt > . 



(c) For each r < A and e > 6 there exists T > such that 

max |6(0) - 6/(0)| < r => max I6(*) - 6(*)l < £ Vt>T. 

i, 3 i, 3 



7 Conclusion 

For nonlinearly coupled harmonic oscillators we have shown that synchroniza- 
tion manifold is semiglobally practically asymptotically stable in the frequency 
of oscillations. Our assumption on each coupling function is that it is locally 
Lipschitz and, if nonzero, its graph lies in the first and third quadrants and 
does not get arbitrarily close to the horizontal axis when far from the origin. 
Our assumption on the interconnection graph is the minimum; that is, it is 
connected. 

One last remark we want to make is the following. If we look at we 
realize that H = [0 1] is not necessary for the rest of the analysis. In fact any 
nonzero H G R lx2 is no worse than [0 1]. For instance, for H = [1 0] coupled 
harmonic oscillators would be represented by 

q t = upi+^jijiqj -qi) 

pi = -uqi 
for which Theorem [3] is valid. 
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